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Abstract. Traditional routes to Charge-Density- Wave in transition metal dichalco- 
genides, relying on Fermi surface nesting or Jahn- Teller instabilities have recently been 
brought into question. While this calls for exploration of alternative views, paucity of 
theoretical guidance sustains lively controversy on the origin of, and interplay between 
CDW and superconductive orders in transition metal dichalcogenides. Here, we ex- 
plore a preformed excitonic liquid route, heavily supplemented by modern correlated 
electronic structure calculations, to an excitonic-CDW order in lT-TiSe2. We show 
that orbital-selective dynamical localisation arising from preformed excitonic liquid 
correlations is somewhat reminiscent to states proposed for d- and /-band quantum 
criticality at the border of magnetism. Excellent quantitative explication of a wide 
, range of spectral and transport responses in both normal and CDW phases provides 

1 strong support for our scenario, and suggests that soft excitonic liquid fluctuations 

mediate superconductivity in a broad class of transition metal dichalcogenides on the 
(N " border of CDW. This brings the transition metal dichalcogenides closer to the bad ac- 

tors in d- and /-band systems, where anomalously soft fluctuations of electronic origin 
are believed to mediate unconventional superconductivity on the border of magnetism. 
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1. Introduction 

Twenty six years after the path-breaking discovery of high-T c superconductivity in 
doped, quasi-two dimensional copper oxides (cuprates), the list of ill-understood strongly 
correlated electronic systems (SCES) with partially filled d or / bands continues to 
grow [H, El E]. The paradigm shift engendered by cuprates also fuelled renewed 
interest and study of older systems. Before the cuprate revolution, these were sore 
spots in the standard model of electrons in metals, centred around the celebrated 
Landau Fermi liquid (LFL) theory. Intense activity to unravel increasing variety 
of unconventional ordered states with ill-understood metallicity has stabilised this 
paradigm shift. Cuprates are not, it is now clear, an isolated example: rare-earth systems 
close to magnetic instabilities at T = 0, increasing number of <i-band perovskites and 
the recent explosion in Fe-arsenides are but a few examples of a truly diverse zoo of 
strange systems. 

Perhaps equally remarkable is the fact that careful work in the recent past has 
brought out unexpected similarities between the newer and older bad actors above. 
This is best exemplified by the recent revival of interest in quasi-2D transition metal 
dichalcogenides (TMD). For more than forty years, the origin of charge density wave 
(CDW) and superconductive (SC) orders and their interplay in the layered TMD had 
remained ill-understood issues jU E]. Historically, appealing to one-electron (density- 
functional theory) band structure led to the conventional wisdom of CDW arising 
from Fermi surface (FS) nesting, or via a band Jahn- Teller (JT) instability. Recent 
revival in the field was stimulated, among other things, by falsification of these views by 
recent ARPES work [6]. The exciting possibility of an alternative, intrinsically strong 
coupling view, involving CDW order emerging as a Bose-Einstein condensation (BEC) 
of an incoherent preformed excitonic liquid (PEL) normal state, arose [7J EJ [9] as an 
attempt to address this new conflict, and has the potential to bring TMD into the 
list of strange actors. Indeed, parallels between TMD and cuprates have increasingly 
been claimed in certain studies [10]. It is also suggested that the absence of magnetism 
in TMD allows studies of a bad metal without attendant magnetic fluctuations that 
complicate the physics in cuprates [TTJ . That a strong coupling scenario for TMD closer 
to the Mott limit is in order, is shown by the fact that the related system lT-TaS2 
has long been understood as a Mott insulator |12j : it is then fully conceivable that 
Mottness also plays important roles in other TMD. From this viewpoint, much as in 
the more recent examples of intense interest, melting of the excitonic-CDW order in 
a PEL scenario would necessarily enhance excitonic fluctuations, which could act as 
an unconventional electronic 'pairing glue' for emergence of the competing SC order. 
However, the theoretical situation still remains unclear, and FS-nesting as well as strong 
electron-phonon views are still claimed to be the instigators of CDW order and SC in 
pressurised or doped systems [TO]. 

The long-studied lT-TiSe2 is a particularly relevant case in point. Careful 
experimental work has unravelled behavior that consistently fails to fit into any 
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conventional views. We list the problems here: 

(i) While optical studies [13] show that the normal-CDW transition is a semi-metal 
to semiconductor transition a la the Overhauser scenario, noticeable spectral weight 
transfer (SWT) over an energy scale of order 1.0 eV upon heating from 10K to 300K, 
along with a clear isosbestic point as a function of temperature (T) below T c d w = 200 K 
and a large normal state scattering rate support sizable electronic correlations. 

(ii) Pure lT-TiSe2 shows bad-metallic resistivity, much above the Mott limit, even at 
pretty low T (<< T c d w ), with a maximum, but no anomaly, in dp/dT at T ~ T c d w [14J. 
Tellingly, recovery of good metallicity at high pressures also destroys SC. Such a 
correlation is repeatedly seen in many systems where unconventional SC appears 
near Mott insulators and magnetism [15] (however, TMD never show magnetism, nor 
has Mottness been hitherto considered important, except for lT-TaS2 p2] and 2H- 
TaSe2p], HI])- Thus, whatever causes strong normal state scattering also facilitates SC 
pair formation below T c . 

(iii) ARPES data [16] shows that the renormalised electronic structure in the normal 
state (T > T cdw ) already resembles an LDA dispersion modified by excitonic correlations, 
supporting a PEL view. Moreover, lack of clear polaron effects (e.g, kinks in 
the band dispersion at phonon energies, well below E^ in resonant inelastic X-ray 
scatteing [8]) and, more importantly, the insensitivity of ARPES band dispersions to 
atomic displacements occuring across T c d w , argue against a band-JT instability, even 
though electron-lattice coupling is ubiquitous to TMD [T7j . 

(iv) Finally, in spite of the incoherent metal features pointed out above, the effective 
mass is only weakly renormalised above T c d w . 

These observations put strong constraints on an acceptable theory. Most 
importantly, they conflict with both FS nesting and band- Jahn- Teller (JT) views 
on general grounds: Absence of band quasiparticles in the normal state rules out a 
conventional ordering instability involving the LDA Fermi surface (FS) features. Simply 
put, the very concept of a well-defined FS in the LDA sense becomes untenable in bad 
metals. As pointed out above, ARPES data also argue against a band-JT instability to 
CDW order. Weak mass renormalisation above T C( j, w conflicts with one-band modeling, 
but not with a multi-band approach jT8]: even in the classic Mott case of V2O3, the 
effective mass in the correlated metal is only moderately enhanced near the (undoubtedly 
correlation driven) metal-insulator transition [THJ EH]- Thus, taken together, these 
observations force one to view emergence of the CDW as a strong coupling Overhauser 
instability of a multi-band bad-metal normal state without LFL quasiparticles. 

These issues have been selectively addressed earlier within conflicting theoretical 
views, but the above observations force us to critically re-examine them. Here, we show 
that these anomalous responses are naturally understood within our PEL view [7] for 
lT-TiSe2, strongly supporting the PEL as a novel (and perhaps generic) alternative 
to conventional theories for TMD. In remarkable parallels with d and /-band systems 
close to Mottness and unconventional order (s), we unravel a new instance of a selective 
metallic state in the PEL, and construct a scenario for the instability of this PEL to a 
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Figure 1. (Color Online) Non interacting DOS, tight binding band structure for the 
Ti-d t2 g and Se-p bands and the corresponding Fermi surface plot. The MO DMFT 
involves Ti-d (black) and Se-p (yellow) bands. 



low-T CDW state. Our philosophy is opposite in spirit to weak-coupling Fermi liquid 
ones, and is actually closer to resonating valence bond (RVB) ideas [20], in the sense 
that ordered state(s) arise as two-particle (BEC) instabilities of an incoherent liquid of 
preformed excitons. 

2. Method 

LCAO band structure for lT-TiSe 2 was constructed [9j [21] by using the Ti-d and Se-p 
states. This gives two bands closest to E F (predominantly Ti-d xy and Se-p z ) as well 
as the FS in very good accord with earlier results [21], as shown in Fig. [U A sizable 
d xy — Pz mixing hybridises the small number of electrons and holes. In this situation, 
even moderate electronic correlations (< 1.0 eV) facilitate exciton formation already at 
high T, as already discussed by Halperin et al. [5] in the sixties. Given that electronically 
active states comprise d xy and p z band states, an intraband Hubbard U ~ 1.0 eV and 
interband correlation U a b — 0.5 — 0.7 eV are realistic values: these can be estimated 
from a first-principles, constrained-LDA calculation, and we believe [22] that these fall 
in the range quoted above. 

Though static Hartree-Fock (HF) theory can now give a CDW ground state [16j [9] , 
description of observed normal state incoherence, and in general, preformed liquid-like 
electronic states, lies outside its scope. Thus, the character and ordering instabilities of 
such bad metals cannot, by construction, be rationalised by appealing to static-mean- 
field theory: this can only be reliably accessed by approaches which can adequately 
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capture dynamical correlations. For lT-TiSe2, a minimal two-band Hubbard model 
as denned below is mandated by LCAO results, and adequate treatment of dynamical 
correlations underlying incoherent behavior is achieved by dynamical-mean field theory 
(DMFT). DMFT and cluster-DMFT approaches have, by now, a proven record of 
successfully treating strong dynamical fluctuatons in correlated electronic systems, 
making them preferred tools of choice in the present context. The two-band Hubbard 
model we use is (calling d xy = a and p z = b) 

Hel = ^2 + ^lm)cl la C kma + U ^ n m n Ui + U ab ^ n ^ n ib 

k,i,m,(T i,l=a,b i 

where l,m run over both band indices a, b, the intra-orbital correlation is U (taken 
to be same for a, b bands, we have checked that results are insensitive to this choice 
within reasonable limits), U a b is the inter-orbital correlation term that, along with 
will play a major role throughout. Further, in TMD, the most relevant A lg phonon 
mode couples to the inter-band excitons [9J by symmetry, and the electron-phonon 
coupling is H el _i = gY,i( A i + A l)(4a c ib + h.c). To solve H = H el + E ei _i within 
DMFT, we have combined the multi-orbital iterated perturbation theory (MOIPT) 
for H d p] with the DMFT for polarons by Ciuchi et al. [231 EI] (see Appendix). 
Actually, this involves extending the polaron-DMFT to multiband cases, seen by 
writing H ei ^i = gJ2A A i + A h( n i,+ ~ n i-) b Y a rotation, c i>± = (c iA ± c i}b )/y/2 and 
rii t ± = c|± c i,±- Finally, as a theoretical advance, we extend the normal state DMFT to 
the broken-symmetry CDW phase as before [7]. This is justified, since, from the above 
discussion and LCAO+DMFT results below, we find, a posteriori, an incoherent PEL. 
Instability to CDW order then cannot occur via the traditional band-folding of well- 
defined Fermi liquid (FL) quasiparticles. Rather, as coherent one-particle inter-band 
mixing is inoperative, ordered states must now arise directly as two-particle instabilities 
of the bad metal. For lT-TiSe2, the residual two-particle interaction, obtained to second 
order is proportional to t 2 ab , more relevant than the (incoherent) one-electron mixing, t a f,. 
The interaction H res ~ -t 2 ab Xab(0, 0) J2<i,j>,aa> c laa c jbac] ba ,c iaa , with Xab the inter-orbital 
susceptibility. Decoupling this intersite interaction in a generalised HF sense yields two 
competing instabilities: H^ F) = - J2<i,j > ,aa>( A ib4 ar7 c iaa + A ab c| a(T c J t fe _ (T + a ->• b), with 
A c dw = (Ai a — Aib) oc (n a — rib) representing a CDW and A a f, cx (ci aa Cjb- a ) a multiband 
spin-singlet SC. Following earlier procedure [7], we compute DMFT spectral functions 
and transport properties in the CDW state at low T, leaving SC for future. 

3. Results 

We now show how the approach envisaged above gives a very good account of a whole 
range of physical responses. We show the DMFT many-body DOS for the 'a' (Ti-d) 
and 'b' (Se-p) bands in Fig. Eh and the corresponding self energies in Fig. Eb, Et bands 
at high (T > T cdw ) and low (T < T cdw ) temperatures. 
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Figure 2. (Color Online) DMFT temperature evolution for (a) DOS of Ti-d and Se-p 
band (ones with peak below FL are for Se-p), (b) Imaginary part of self energy and 
(c) real part of self energy for Ti-d band and Se-p band. The self energies with a dip 
like structure close to -2 eV are for Se-p band. 
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The DOS shows clear 'semi-metal' features, which, however, are not those of a 
conventional semi-metal. Examination of the self-energies bares the hidden selective 
Mottness in the system: Im£ a (w,T = 200 K) clearly shows incoherent metal features 
(ImE (a> = 0) > at E F {= 0)), in accord with (dp/dT) < for high-T (black dashed 
curve in Fig.2b), while the 6-band continues to show renormalised band insulator-like 
features at low energy. This orbital-selectivity (OS) is thus related to strong normal 
state scattering, rather than the onset of CDW order. Finding of OS in lT-TiSe2 is 
surprising (Ti being nominally d°), but can be traced back to the fact that, in presence 
of a small number of electrons and holes induced by t ab , even a moderate U a b leads 
to exciton formation, and thence to selective-Mottness in the multi-band situation 
that obtains in lT-TiSe2- Recall that the original idea of Mott was in fact indelibly 
tied to quantum melting of correlation-induced excitons under pressure. Intriguingly, 
in lT-TiSe2, stabilisation of excitonic CDW order enhances selective- Mott features at 
lower T, as emergence of the co = pole in ImS a b (u;) clearly shows. Simultaneously, 
however, examination of ReE Qj b(u;) also clearly shows small mass enhancements. In 
extant literature [8], this finding below T cdw is proposed as a concrete example of a 
reduced effective mass due to renormalisation effects caused by the appearance of a 
strong periodic potential below T c ^. Our results are certainly in accord with this 
finding. As mentioned before, however, no large mass enhancements are theoretically 
(at least within DMFT) expected in a multiband system, even in an orbital-selective 
Mott regime. This goes hand-in-hand with slight stabilisation of the normal state gap, 
now interpreted as a true CDW gap. As advertised before, and in full accord with data, 
sizable T-driven SWT accompanies the transition. 

If the PEL alternative is to be credible, the full range of observations must be 
explicable without any additional assumptions. In Fig. [3] we compare our DMFT 
one-particle spectral function, A a) &(k, u) = — ImG ajb (k, uj)/tt, and renormalised band 
dispersion, -Ek,a,b = ek,a,b+ReSk, a ,b, with ARPES data [IE] where these exist. This 
is already a stringent test for theory: while LDA plus static-HF can conceivably 
yield agreement with band dispersions, the real test is a simultaneous description of 
the ARPES lineshapes. For a dynamically fluctuating 'liquid', the latter is expected 
to show broad continuum-like features without infra-red (LFL quasiparticle) poles. 
Rather remarkably, our LCAO+DMFT results provide an excellent semi-quantitative 
description of extant ARPES dispersions and lineshapes up to rather high energy. In 
particular, they bare the preformed excitonic features in i?k,a,b and clear, associated 'gap' 
features in ARPES lineshapes above T cdw . Excellent accord in all details, including the 
band positions, their intensity distributions, and band-shifts as a function of T, with 
the ARPES dispersion is clear from a direct comparison between our Fig. H^a-c) with 
Fig. 2 of Monney et al. [16] In particular, we can even identify the detailed features in 
the T-dependence of the ARPES intensity (Fig. [3]) with data: {%) the peak at ~ —0.2 eV 
is identified as the valence band backfolded to the M point by comparing the red curve 
in Fig. [3] with the top panel in Fig. HI (ii) a new peak, labelled 'C by Monney et 
al, also appears below T cdw and gets more pronounced upon decreasing T, precisely as 
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Figure 3. (Color Online) Theoretical ARPES at M point at different temperatures 
(T). As T increases across T c d w , the peak in the DMFT spectral function crosses the 
renormalised Fermi energy (Ep). 

seen. (Hi) The peak '£)', identified as a 'quasiparticle peak' originating from coupling to 
phonons in experiment is also obtained in DMFT in very good qualitative accord with 
data: especially interesting is that it is quickly damped out as T increases, exactly as 
seen in the top four curves in our Fig. [3l However, at high T, we resolve an additional 
(sharper) low-energy peak, also arising from electron-phonon coupling, which remains 
dispersionless above T C( i w : this feature is not seen in the ARPES results of Monney et 
al. Finally, we also cannot observe peak i A\ identified with a second, spin-orbit split, 
valence band, since we have not used the full set of Se-p and Ti-d bands in this work. 
Nevertheless, the accord between theory and experiment is excellent with regard to 
features relevant for the CDW in lT-TiSe2- 

We now consider the ARPES 'band structure' maps. At high-T, along K-M-K 
direction in the Brillouin zone the putative valence band, (VB) heavily damped by strong 
(due to preformed excitons) scattering, crosses E F {= 0) while the putative conduction 
band (CB) lies totally above Ep. This conflicts with the LDA (or LCAO, see Fig. [T]) 
results, which predict a sizable overlap between VB and CB states at Ep, and reflects the 
failure of the FS nesting mechanism in the real one-particle spectra (in fact, to cure such 
conflicts, the standard explanation within conventional views has been to attribute such 
discrepancies between LDA and ARPES data to uncontrolled excess [8] of Ti. Our view 
naturally produces this, without having to take recourse to such extraneous conditions). 
Clear Hubbard-band shake-up features in DMFT are also seen in the color-plot: in 
particular, we predict that the u > part of the lineshape (at M point in Fig. 0H) should 
lend itself to observation in inverse ARPES (ARIPES) studies in future. In addition, 
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Figure 4. (Color Online) Theoretical ARPES map along K-M-K direction at (a) 20K 
(b) 150K (c) 270K (d) Total EDC at M point at different temperatures. Excellent 
accord with ARPES data is clearly seen. 
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comparison of the DMFT energy distribution curve (EDC) with extant ARPES result 
at the M point also shows (Fig.3) very good accord as a function of T over the whole 
range. A moderate band narrowing of LCAO bands and broad spectral lineshapes with 
appreciable T-dependent SWT in one single system are generic fingerprints of strong 
dynamical correlations in multi-orbital systems, and so, as alluded to before, the second 
feature above, essential for describing 'liquid' correlations, cannot be accounted for by 
a static-HF theory, as done so far j9]. 

Detailed Fermi surface (FS) maps as a function of T in lT-TiSe2 are rare. We have 
taken the FS mapped out by Rossnagel et al. [25] to study how the PEL idea survives 
this important test. While instabilities driven by LDA FS forms the backbone of weak 
coupling or itinerant views, it is also possible, in a renormalised itinerant or Mottness 
based theories, that these could involve a new FS sizably reconstructed by correlations. 
That the FS does not reconstruct across T c d w in lT-TiSe2 was one of the main (among 
others, see above) arguments for invoking unconventional PEL scenarios in the first 
place. Close inspection of FS evolution across Td W shows that, while this is undoubtedly 
correct, there are still specific features which any theory needs to confront: (i) the band 
pockets are smeared out at high temperature (T > T c d w ) and, more importantly, a much 
brighter ring structure appears at the M point below T cdw . 

In Fig. [51 we show our DMFT FS at low (top panel) and high-T (bottom panel). 
As expected due to strong PEL scattering, the high-T FS is sizably smeared out (arises 
from the finite Im£ a (a; = 0) at 200K from Fig. |2J). Remarkably enough, we also find a 
clear ring feature at the M point below 150K (lower panels), in remarkable accord with 
data. Since this ring-like feature becomes well-defined only below T cdw , it is a direct 
consequence of CDW order-induced reconstruction of electronic states. This is fully 
consistent with the excitonic CDW view, where CDW follows excitonic 'solid' order, 
and the latter arises principally via interplay between t^ b and U a b, implying backfolding 
of Se p band from the T point due to the 2x2x2 CDW superstructure formation. Further, 
close observation of DMFT results shows that the central pocket around the T point 
has seemingly developed 'elongated' shape instead of the hexagonal shape expected 
from the LCAO result. We believe that this important modification, hitherto not noted 
sufficiently, arises due to an orbital- dependent electronic structure reconstruction that 
is essentially driven by the k-dependent form factor of the inter-orbital hybridisation 
(t a b(k)). Experimental confirmation of this feature would thus constitute additional 
support for a PEL scenario. Finally, the stabilisation and slight increase of the normal 
state 'gap' below T c ^ also accords with optical data [13] and with the semi-metal-to- 
semiconductor characterisation of the normal-CDW transition in lT-TiSe2. 

Thus, such quantitative agreement between our excitonic-DMFT results premised 
upon a novel PEL view and ARPES in all important details lends strong credence to 
the idea of a dynamically fluctuating excitonic liquid at high T giving way to a low-T 
CDW order. However, to further qualify as a credible candidate, the same formulation 
must also describe transport as well. Fortunately, in DMFT, this task is simplified: it is 
an excellent approximation to compute transport co-efficients directly from the DMFT 
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Figure 5. (Color Online) DMFT FS map at (i) 20K (ii) 150K (iii) 270K, showing very 
good agreement with the T-evolution of the Fermi surface in ARPES data of Rossnagel 
et al. [25] 
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Figure 6. (Color Online) DMFT results for T-dependent optical conductivity, showing 
good agreement with optical data across T c d w . 

propagators G at b(k,u) [26J, since (irreducible) vertex corrections rigorously vanish for 
one-band models, and turn out to be surprisingly small even for the multi-band case. 

In Fig. |6l we show the optical conductivity cr(u) as a function of T, wherein very 
good accord with data up to an energy ~ 0(0.8) eV is clear. A clean CDW gap 
at low T closes rapidly with increasing T via rapid spectral weight transfer from the 
relatively sharply defined hump at 0.4 eV to the infra-red regions, precisely as seen. 
Given that we have kept only the two bands crossing Ep in LCAO, agreement at higher 
energy (> 0.8 eV) is not expected. The relative sharpness of the low-energy optical 
response at low T is deceptive: it is not a FL Drude peak, and in fact, (fully consistent 
with the selective-Mott behavior in the DMFT spectra above) no FL coherence sets 
in, even at lowest (10 K) T in DMFT. It is rather a reflection of reduced incoherence 
due to CDW gap opening. Finally, we also resolve an isosbectic point in a(u,T) curves 
(where a(0.2eV) remains invariant) at different T: this is clear manifestation of sizable 
dynamical correlations [7], and is again in good accord with the isosbectic point seen 
around 0.27 eV in the optical study [13]. 

The DMFT resistivity (Fig. [7]) also shows an insulator-like behavior above T C( } W , 
a broad peak without any anomaly below T c d w = 200 K, and bad metallic behavior 
below Td w , in full accord with data [14]. The insulator-like behavior for T > T cl i w is 
a reflecting of the large normal state incoherence, whence carriers strongly scatter off 
fluctuating, incoherent (preformed) excitons. The bad-metallic behavior below is 
then naturally attributable to reduction of this strong normal state scatttering due to (i) 
opening up of a CDW gap, and (ii) concomitant increase in the tendency to excitonic 
coherence. Given the bad-metallic resistivity, the consequent short scattering mean- 
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Figure 7. T-dependent dc resistivity within DMFT. Except for the low-T part, 
LCAO+DMFT results agree very well with experimental results, including a broad 
maximum, rather than any sharp non-analytic change, across T cc i w (see text). 

free path (in fact, hpl ~ 0(1)) invalidates a quasiclassical Boltzmann equation-based 
approach to transport. Interestingly, this is precisely the regime where DMFT should 
work best. The absence of a sharp ordering anomaly at T c & w is additional evidence 
against a weak coupling view of the instability. In fact, in a weak-coupling instability, 
resistivity should have shown a sharp ordering anomaly at T c d w on general grounds [27] 
(in addition, dp/dT must also show critical behavior, with exponents linked to those 
extracted from thermodynamic measurements). The strong coupling view is further 
supported by finding of a large 2A//c B T C(Ju , ~ 7— 10 in TMD (about 7 for lT-TiSe 2 and 
10 for 2H-TaSe2 [21]). This feature is reminiscent of high-T c cuprates [29] and implies 
that CDW formation is associated with a BEC, rather than a BCS-like scenario for the 
exciton instability. Also, strong inelastic scattering needed to rationalise bad-metallicity 
above T c d w and sizable T-induced SWT are characteristic signatures of a strong coupling 
limit. 

Thus, taken together, very good accord with ARPES and transport data strongly 
supports our basic hypothesis: {%) the normal state is a strongly fluctuating liquid of 
incoherent excitons, and (ii) CDW order in lT-TiSe2 must fall into the strong coupling 
class, qualitatively different from a conventional Overhauser transition of well-defined 
band (LFL) quasiparticles. 

Emergence of CDW order from a PEL should leave further specific signatures in 
other data. First, CDW and excitonic correlations must now track each other beyond 
T c dw j wen hrto the PEL state. In Fig. [Sj, we show the excitonic and CDW order 
parameters, A exc oc (T x ) = ({c\ a c ib + h.c)),A cdw cx (T z ) = (n a - n b )/2, along with their 
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Figure 8. (a)Expectation values and derivative (inset) of inter-orbital excitonic 
(ImGab), and of (b) CDW order parameter, n a -nb = A cc [ w and its derivative (inset) as 
a function of T. 



T-derivatives, as a function of temperature (T). In particular, as expected in the PEL 
scenario, both follow each other and are finite way above T c d w , in nice qualitative accord 
with the T-dependence of the CDW order parameter from ARPES [16] . The absence of 
a conventional (BCS-like) mean-field transition can be rationalised by noticing that, in 
presence of the el-lattice coupling, H e i-i = gJ2i(Ai + A\)(cl a cn, + h.c), the universality 
class of the normal-CDW transition turns out to be that of an Ising model in an external 
field. This implies that the (mean-field) transition is smeared into a smooth crossover. 
The relevance of excitonic correlations is also visible as follows. We see that p(T) 
closely follows the T-dependence of dA exc /dT above (Fig. E]), but does not show any 
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Figure 9. (Color Online) Phonon spectral function at different temperatures. The 
inset shows change in peak position of the A\ g phonon mode with T, in good accord 
with Raman data above T cc i w (see text). 

obvious correlation with dA c d w /dT: this shows that the anomalous resistivity is caused 
by carriers scattering off incoherent excitonic (liquid-like) correlations. In our PEL 
scenario, the change to metallic behavior in p(T) is attributable to reduction of the 
strong scattering when excitons condense at T c ^ w . True CDW order now follows a true 
BEC of preformed excitons at lower T. 

Given the specific form of the electron-phonon coupling, now also interpretable as 
coupling of interband excitons to A lg phonon mode by symmetry, the lattice is expected, 
very generally, to react to the T-dependent changes (incoherent at high T, more coherent 
at lower T) in exciton dynamics. In particular, the DMFT phonon spectral function, 
computed from p p h(uj) = (—l/n)lmD(u) = (— l/7r)im[2fi /(u; 2 — Q% — 2g 2 VLoXab{^))} 
and shown in inset of Fig. [9j should mirror excitonic CDW correlations. 

Several interesting features stand out: (i) p p hiy>) shows maximum intensity with 
reduced linewidth around T m = 150 K, precisely where the resistivity peaks, (ii) Above 
T cdw , the phonon spectrum is noticeably broader, reflecting coupling to incoherent 
excitonic liquid modes, and (iii) asymetry in p p h{u) increases as T is lowered below 
T c dw This reflects precursors of a Fano-like structure, arising from treatment of H e i~i 
beyond the adiabatic limit within our DMFT. Onset of excitonic coherence below T c d w 
reduces strong normal state scattering, sharpening the phonon spectrum, (iv) Finally, 
we find that the detailed T-dependence of both the y4 l5 -phonon lineshape and linewidth 
excellently matches that extracted from a Raman scattering study [30]. Only below 
T c d w is there a discord between measured and computed linewidths: it levels off to a 
constant in experiment, a feature not captured by our present calculation. The reason 
is that we have not carried out a study of the change in phonon dynamics below T c a w . 
A full study of lattice effects must rest on a more realistic input to the full phonon 
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spectrum of lT-TiSe 2 , a point we leave for future study. Thus, taken together, very- 
good quantitative accord with a whole host of spectral and transport responses for 1T- 
TiSe2 and a comprehensive qualitative rationalisation of structural features in one single 
theoretical picture constitutes overwhelming support for a novel PEL view. Thus, the 
central conclusion of our work is: The CDW instability must now be interpreted as a 
strong coupling instability of an incoherent PEL normal state, rather than as a weak 
coupling Fermi surface nesting instability of a good LFL. 

The above findings have important implications for SC arising under pressure. 
A finite pressure increases a — b band overlap via t a b, leading to a redistribution of 
electrons and holes, and weakens excitonic CDW order. From Fig. this implies 
enhancement of excitonic fluctuations, also seen by noticing that a fall of (r|T|) must 
transfer weight to the transverse part, i.e, enhance ((r^Tj + h.c)) via the pseudospin 
sum-rule T 2 = T(T + 1). Thus, near the critical p = p c where CDW is destroyed, 
we expect a maximisation of excitonic fluctuations. These can constitute the critical 
collective (electronic) fluctuations that could lead to an instability to an inter-band 
SC (competing with CDW) with a finite A ab ~ (ciavCjb-a) ■ However, normal state 
incoherence implies that SC must now arise directly from the incoherent metal via a 
two-particle instability: this can be explored by solving H = H e i + H res in the pair 
channel. The el-lattice coupling will further enhance the effective pair interaction in 
H res above [7J. Even without the benefit of a detailed DMFT calculation for the SC 
phase, we thus expect the SC T c to be maximal close to p c if dynamical excitonic 
fluctuations mediate cooper pairing: this is again in accord with observations p3] and 
the causal link is compelling. 

We believe the present work also has further important implcations on a broader 
level. Our unexpected finding of OS bad-metallicity in the high-T PEL brings 1T- 
TiSe 2 closer to the bad actors in d- and /-band systems, where OS is a consequence 
of orbital-selective Hubbard correlations [31] (or, in the cuprates, momentum-selective 
correlations [29]). Given that lT-TaS2 has long been known to be a Mott insulator, 
our proposal of selective-Mottness and associated excitonic liquid features is quite 
reasonable, albeit novel. However, since magnetism is never an issue in the TMDs 
of interest, what underpins anomalous responses reminiscen of d- and /-band critical 
systems near magnetism is moot. The common unifying mechanism seems to be 
selective-Mottness, which is indeed one of the main contenders for understanding 
anomalous quantum criticality in d and /-band systems [321 E3]- Selective-Mottness 
itself can always occur via different mechanisms: in particular, depending upon the 
microscopies of a given system, it can be associated with magnetism, but it clearly does 
not always need to be so. It can also be associated with other density-wave instabilities 
in the charge or orbital sectors in mutliband systems such as TMDs. Nevertheless, 
a common element, namely destruction of Landau FL quasiparticles, would always 
accompany onset of selective-Mottness. Thus, in light of the present work, it should not 
be too surprising that the high-T PEL bears some similarities with the FL* theory [32] 
in the /-electron QCP context. Evidence attesting to this comes from the fact that 



Preformed Excitons, Orbital Selectivity, and Charge-Density-Wave Order in !T-TiSe2l7 

only the Se-p band crosses Ep, while the Ti-d band lies above Ep already in the 
high T liquid. In other words, the observed FS has already reconstructed from its 
LDA counterpart to reflect the preformed excitonic character of the normal state. 
Spectral and transport features show characteristic incoherence features, as expected 
in an FL* state, and the efficacy of DMFT as a valuable tool to understand these is also 
known in the /-electron context [331 ET] - That this identification implies an interesting 
scenario, similar to those invoked for anomalous metals on the border of (T = 0) 
magnetism may also permit a rationalisation of the apparent 'similarities' claimed 
to exist between under-doped cuprates and TMD [10]. Viewed from the perspective 
of our study, these similarities are, ultimately, manifestations of the selective- (bad) 
metallicity in both cases: momentum-selective in under-doped cuprates, and orbital- 
selective in lT-TiSe2. Thus, we arrive at a perhaps general idea with broader appeal 
across classes of correlated systems: Selective Mottness (whether orbital or momentum) 
induces critical electronic liquid states characterised by loss of LFL coherence. This 
collectively fluctuating electronic fluid can subsequently become unstable, either to a 
myriad (charge, spin, orbital) of density wave orders, or to (competing) unconventional 
superconductive orders, depending upon the actual microscopies of the system under 
study. In TMD, material specific reasons favor competing inter-band CDW and SC 
orders from such a strongly fluctuating excitonic liquid. 

4. Conclusion 

In conclusion, we show that a whole host of physical responses in lT-TiSe2, difficult to 
reconcile with a band FS nesting or band-JT mechanisms, is provided a natural and 
good quantitative explication within a new PEL alternative. Finding of novel, orbital- 
selective Mott features in renormalised normal state electronic structure via DMFT 
leads to an excellent description of a range of spectral and transport data in both nor- 
mal and CDW states. More importantly, it brings the TMD closer to the anomalous 
critical metals of much more recent interest. Along with its success for 2H-TaSe 2 [TJ, 
our work relates the PEL idea to a more generic theoretical level for TMD. Particularly 
interesting should be to test how this new proposal fares for the even more correlated 
lT-TaS 2 [35], which is undoubtedly a known Mott insulator [12], as well as competition 
between CDW and SC, in future work. 
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helpful discussions. 
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5. Appendix 

5.1. Tight binding band structure 

TiSe 2 is a layered material with hexagonal layers of Ti sandwiched between layers of Se 
atoms: a lT-polytype transition metal dichalcogenide. The Ti d-band is in a d° state 
in TiSe2- Due to octahedral coordination the Ti atom d-shell is split into a set of high 
energy e g orbitals and low lying, degenerate t2 9 orbitals. We consider only the low lying 
t2 9 orbitals, for, they are the ones that form the hybridized bands with Se-p orbitals 
close to the Fermi level. In the LCAO calculation, therefore, we consider charge transfer 
between them and the surrounding Se 4p orbitals, so that six Se p orbitals and three 
Ti d orbitals per unit cell are involved. Our results are similar to those of Wezel et 
al. [9], but we do not invoke quasi-lD features resulting from the orbital dependence of 
the hopping matrix elements as the driving cause of CDW. However, the shape of the 
central pocket around the T point in the DMFT Fermi surface in the main text does hint 
toward such anisotropic features getting more pronounced in the CDW ordered phase. 
The tight binding Hamiltonian is thus constructed as 

H = ^ K\ A« - Pli,aPli,a ~ Pli,J 2 i,a) 

<i,j>,a,/3 

<i,j>,a,/3 

t^Mi^W + H.c.}). (5.1.1) 

Here d^p^ and p' 2i create electrons on the Ti, the upper Se and the lower Se atom 
respectively. The labels a and run over all possible orientations of the Ti d-t 2g and 
Se p orbitals and < i, j > denotes neighbouring sites. t dd ,t pp and t pd are the different 
hopping matrices for d and p orbitals and A is the chemical potential. 

This Hamiltonian results in a 9x9 matrix using Slater-Koster integrals, and several 
matrix elements vanish rigorously due to crystal symmetry. We obtain the LCAO band 
structure by diagonalizing the resultant matrix. By adjusting the values of Slater-Koster 
integrals and the chemical potential the tight binding bands can be fit to the extant 
LDA calculation. Setting dda = — 0.2, ddir = 0.2, dd5 = —0.5,pd7r = 0.5, ppa = 0.4 and 
A = 2.0 we get a TB fit in good quantitative agreement with the LDA results J2Tj [55] . 
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5.2. Electron phonon Self Energy 

To calculate electron-phonon self-energy we incorporate the procedure first used by 
Ciuchi et al. [23] (using Einstein phonons) into our multi-orbital DMFT. The local 
intra- and inter-orbital Coulomb correlations are treated upto second order selfconsistent 
multiband IPT as usual. Since both Hubbard and electron-phonon interaction terms are 
local, their combined effect can be treated simultaneously within DMFT. Contribution 
of the electron-phonon coupling to electronic self energy is given by 

f E GMDM = A "'I nAC i + + 1 ~ "'M ] (5.2.1) 

ILU„ 

Where N q = 

The full Hamiltonian, including excitonic coupling of phonons is 

H= ^2 {tt+e a 5 ab )cl ar7 c kba +U ^ n i^ni^- a +U ab ^ni a ni b +g^(cl^c ba +h.c.)(Al+Ai 

k,a,b,a i,fi=a,b i i 



-y E c ^(i - +u J2 4 a. (5.2.2) 



(2) 

We start with an initial ansatz for the self energy T, int (uj) = Un + AT, '(to) where 
Ti^\oo) is the second order contribution of electron-electron and interband excitons 
coupling to Ai g phonons. The IPT scheme is that we calculate lattice Green's function 
(Gfa) from this full self energy, from which the bare Green's function is found via the 
Dyson equation: G^ 1 = Gj^ + E a . Plugging this Go a back into the IPT, we obtain a 

(2) 

new estimate of Sg a and this procedure is iterated to convergence. 

Having both Hubbard-like and el-ph couplings changes the estimate of A ab used 
in the interpolative self-energy in IPT as follows. Following usual procedure, A ab is 
calculated from the condition that it reproduce the leading behavior of the (of the 
exact atomic limit) self-energy at high frequency. The leading behavior for large 
u can be obtained by expanding the Green function in a continuous fraction 



Gf a (k,uj) = l/(u — €f a — M\ a — M2 ^ + Mla ) • Here, denotes the i th order moment of the 
density of states. One can compute these quantities by evaluating a commutator [38J, 
and for the model above, M 2a -M 2 a = U 2 b (n fa (l-2n fa ) + < n fa n fb >)+g 2 (n fa+ +n fa _). 
Here, n/ a is the number density calculated from the full Green's Function, nf a+ = 
9 2 E^( G /( W + u q )(N q + n f (u)) and n /a _= g 2 ^ w (G f (u ~ u q)( N q + n f(~ UJ )). From the 
large frequency limit of (1), Sg a (a;) = U 2 b n 0a (l - n 0a ) + g 2 (n 0a - + n 0a+ ). Here, n Qa is a 
fictitious number density of the 'bare' Green's function. Explicitly, n 0a+ = g 2 ^2 u {Gq(uj-\- 
Uq)(N q + 7if(u))) and noa-= g 2 X] w (Co(w- co q )(N q + nf(— oj)). Comparing with the exact 
high-frequency limit, we thus have A ab = ^ ( "^^:^£^^ +nf "- ) • 

Within the DMFT approximation, the phonon self-energy turns out to be H(q, to) = 
' ; v ,,/ — where Xc(q>w) is the usual charge-charge response function [3H] and 



1+9 Xc(q,u)D (q,u) 



Preformed Excitons, Orbital Selectivity, and Charge-Density- Wave Order in !T-TiSe 2 20 

D (g, uj) is the bare phonon Green's function. Xc(q,u) is estimated by a renormalised 
bubble contribution of the DMFT Green functions. We ignore the irreducible vertex 
corrections, since their contribution should be small for the small number of electrons 
and holes that characterise the two-band system close to an excitonic insulator/liquid 
regime, and is a further approximation. 
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